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ABSTRACT 

In  this  paper  we  study  the  Riemann  Problem  for  a  system  of  conservation 
laws  which  exhibit  internal  friction  similar  to  that  seen  in  viscoelastic 
solids  of  the  maxwell  type.  The  solutions  we  obtain  have  a  single  shock  and  a 
single  contact  discontinuity  and  off  of  these  singular  curves  they  are 
smooth.  The  results  we  obtain  are  two  fold.  First  we  show  this  problem  is 
globally  solvable  in  time;  this  reguires  precise  a-priori  estimates  for  the 
solution  off  of  the  singular  curves.  Secondly,  we  obtain  asymptotic  or  large 
time  information  about  the  solution  which  guarantees  that  in  a  weak  sense  it 
converges  to  special  traveling  wave  solutions  of  the  eguations  with  compatible 
data. 


AMR  (MOS)  Subject  Classifications:  35B40,  35B45,  35L65,  35L67,  73F99 
Key  Words:  Shock  waves,  Riemann  Problems,  viscoelasticity 
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SIGNIFICANCE  AND  EXPLANATION 


Problems  arising  in  continuum  mechanics  are  often  modelled  by  systems  of 
conservation  laws.  For  perfect  materials  such  as  an  elastic  solid  or  an  ideal 
gas,  these  balance  laws,  when  adjoined  with  the  constitutive  equations 
describing  the  material,  lead  to  systems  of  nonlinear  hyperbolic  partial 
differential  equations  in  which  the  characteristic  speeds  are  dependent  on  the 
amplitudes  of  the  motion.  Such  systems  have  the  property  that  nonconstant 
disturbances  are  amplified  and  solutions  which  were  initially  smooth  develop 
discontinuities  in  finite  time.  It  is  well-known  that  this  loss  of  regularity 
can  be  prevented  if  viscous  frictional  forces  are  incorporated  into  the 
constitutive  assumptions  describing  the  material. 

A  different  situation  obtains  when  the  constitutive  assumptions  of  a 
perfect  material  are  modified  to  account  for  long  range  memory  effects  but 
viscous  forces  are  neglected.  The  following  results  are  typical:  (i) 
solutions  generated  by  small  amplitude,  smooth  data  persist  for  all  times  and 
decay  to  a  rest  state  as  time  proceeds  to  plus  infinity,  and  (ii)  solutions 
generated  by  large  amplitude,  smooth  data  develop  discontinuities  in  finite 
time.  In  a  word,  these  results  show  that  for  small  data  the  damping  generated 
by  the  memory  effects  competes  favorably  against  the  destabalizing  mechanism 
generated  by  the  nonlinearity  in  the  system  while  for  large  disturbances  the 
reverse  is  true. 

It  is  also  well-known  that  such  memory  like  continua  support 
nonequilibrium  solutions  called  traveling  waves.  These  solutions  have  a 
richer  structure  than  their  counterparts  for  perfect  materials.  In  the  latter 
case,  such  traveling  waves  are  piecewise  constant  solutions  separated  by  a 
shock  wave  advancing  at  constant  speed.  What  is  not  known  is  whether  these 
traveling  wave  solutions  are  stable  for  memory  like  continua. 


The  purpose  of  this  paper  is  to  study  a  model  system  consisting  of  a 
single  conservation  law  adjoined  to  a  constitutive  equation  of  the  memory 
type.  This  system  supports  traveling  waves  similar  in  structure  to  those  seen 
in  real  systems  of  actual  interest.  What  we  are  able  to  show  is  that  the 
solution  of  this  system  with  piecewise  constant  initial  data  (the  Riemann 
Problem)  converges,  albeit  in  a  weak  sense,  to  the  traveling  waves.  In  the 
process  of  establishing  this  result  a  number  of  techniques  are  developed  which 
should  be  helpful  in  examining  more  realistic  systems. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive- 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report.  ^ 
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1 .  Introduction 

In  this  note  we  consider  the  system 

u*.  +  0  =  0  ,  and  (H) 

t  x 

(a  -  f(u))t  +  (o  -  pf (u) )  =  0  (C) 

where  0  <  p  <  1  and  f  satisfies 

f(0)  =  0,  0  <  f'(u),  and  n  <  f"(u),  0  <_  u  .  (1.1) 

Other  assumptions  will  be  imposed  on  f  as  they  are  required.  This  system  serves  as  a 

model  for  the  equations  of  motion  of  a  viscoelastic  solid.  In  this  analoqy  (B)  represents 

the  balance  law  and  (C)  the  constitutive  eouation ,  In  fact  we  take  (C)  to  be  the  same  as 

the  equation  representinq  a  nonlinear-maxwell  solid  (1,2).  The  relation  of  the  system  (B) 

and  (C)  to  the  equations  of  motion  for  a  maxwell  solid  are  the  same  as  the  relation  of  the 
2 

Hopf  equation,  ut  +  (y- )  -  0,  to  the  equations  of  motion  for  an  ideal  qas;  namely  in 

x 

certain  situations  the  respective  reduced  equations  approximately  describe  what  is 
happeninq  in  the  full  system  and,  as  importantly,  they  provide  one  with  some  intuition 
about  the  more  complicated  systems. 

nur  interest  is  in  the  Riemann  Problem  for  (B)  and  (C);  that  is  in  solutions  of  (n) 
and  (C)  which  assume  the  followinq  data  at  t  =  0: 

f  (u  ,pf (u  ) ),  x  <  0 

(u,0  )  (x,0)  =  <  0.2) 

|  (0,0),  x  >  0  . 

Here,  u  >0.  This  system  is  hyperbolic.  It  has  one  nondeqpnerate  characteristic  field 


's.d.N.Y. -Buffalo  and  the  National  Science  Foundation, 
^Brown  University  and  Academia  sinica  of  china. 


Washinqton, 


P.c. 


Sponsored  in  part  hv  the  united  States  Army  under  Contract  no.  nA'\0?'>-'lP_c- 
nrna  1  ,  and  the  National  s-ienoe  cmm,!,-,.  i  m.  or  lor  r',r**nt  *:o,  ‘m'S-R'  1 1  S 1 1  . 
work  was  completed  while  rreerberq  was  visit.ino  t'-e  ‘lather- itios  nese->roh 
Center  at  the  University  of  sronsin-Madi  sen . 


which  propagates  with  speed  f'(u)  and  one  linearly  deqenerate  field,  namely  the  lines 
x  =  constant.  In  regions  where  (u,0)  is  smooth,  (B)  and  (C)  are  equivalent  to 


o  +  f'tujo^  +  (1  -  u)0  -  p$  «  0,  and  (0) 

♦  t  -  <1  -  +  !*♦  =  0  ,  ($) 
where  t,  0,  and  u  are  related  by 

$  =  f(u)  -a  or  equivalently  u  =  f  '(0  +  $)  .  (u) 

The  system  (B)  and  (C)  also  supports  traveling  waves  which  assume  the  data 
(u  ,  yf(u  ))  and  (0,0)  at  x  equal  minus  and  plus  infinity  respectively.  Since  these 
waves  will  play  an  important  role  in  what  follows  we  record  their  properties  now.  A 
traveling  wave  which  meets  the  aforementioned  boundary  conditions  is  a  solution  of  (B)  and 
(C)  which  is  a  function  of  t  =  t  -  x/c  where 


Uf(u  ) 


(1.3) 


The  u  component  of  this  solution  satisfies 

~  ( f ( u )  -  cu)  =  (cu  -  Pf(u))  (1.4) 

and  the  limit  relations 

lim  u(r)  =  Q  and  lim  u(T)  =  u_  .  (1.5) 

T  -*■-**  T++® 

J  and  $  are  obtained  from  u  by 

0(T)  =  cu ( T )  and  t(T>=f(u(T))-cu(T).  (1.6) 

The  results  for  traveling  waves  are  summarized  below: 

(a)  When  c  <  f'(0)  there  is  a  unique  (to  within  a  translation),  smooth,  strictly 
increasing  solution  of  (1.4)  satisfying  (1.5). 

(b)  when  c  >  f'(0)  there  is  a  nondecreasing  weak  solution  of  (1.4)  which  is 
unique  to  within  a  translation.  The  particular  solution  with  a  jump 
discontinuity  at  r  =  0  has  the  following  additional  properties: 

(i)  u(T)  3  0,  -*  <  T  <  0  , 

(ii)  lim  u ( r )  =  u,  where  u.  >  0  is  the  unique  solution  of 


f(u„)  -  cu#, 


and 


(1.7) 


(ill)  u(t)  satisfies  (1.4)  on  (0,“)  and  meets  the  boundary  condition  (1.5) 
at  plus  infinity. 

The  global  (in  time)  existence  question  for  the  Riemann  Problem  (B),  (C),  and  (1.2)  i3 
easily  resolved.  One  may  either  use  a  fractional  step  variant  of  Glimm's  method  (3)  or  the 
finite  difference  approach  used  by  Greenberg  (4)  for  a  system  of  integro-dif ferential 
equations  of  which  (B)  and  (C)  are  a  special  cast.  Using  an  iteration  scheme  and  the 
contraction  mapping  principle  it  is  also  relatively  easy  to  show  there  is  some  time 
T  ax  >  0  such  that  u,  a,  and  $  have  the  following  properties: 

(i)  (u,or,$)  -  (u_,pf(u_),(1  -  u)f(u_))  for  x  <  0  and  0  _<_  t  <  Tmax  ; 

(ii)  there  is  a  c'tO.T  )  curve  x  »  s(t)  with  xr  >  0  sue h  that 

IndX  Qt 

(u,<J,$)  is  C1  and  satisfies  (B),  (C),  (o),  ($),  and  (u)  on 

0  <  x  <  s(t)  and  0  <  t  <  Tmax, 

(iii)  <u,o,$)  a  (0,0,0),  x  >  s(t)  and  0  _<_  t  <  Tmaxi 

( i v )  the  line  x  =  0,  0  <  t  <  Tmax,  is  a  contact  discontinuity  and  the 

following  limit  relations  obtain: 

lim+  (u,o,$)(x,t)  «  (uQ(t),uf (u_),40(t))  (1.8) 

x*0 

where  uQ(*)  is  defined  Implicitly  by 

f(uQ(t))  -  uf(u_)e"‘lt  +  f(u_)(1  -  e'Ut)  (1.9) 

and 

♦  Q(t)  -  (1  -  u)f(u_)(1  -  e_|Jt)  ;  (1.10) 

(v)  the  curve  x  »  s(t),  0  <_  t  <  Tmax,  is  an  admissible  shock  wave  and 

satisfies  the  Rankine-Huqoniot  conditions  for  (B)  and  (C);  namely  the 
relations 

lim  (u,o,$)(x,t)  -  (U(t),f (U(t) ),0)  ,  (1.11) 

x*s(t) 
x<  s ( t ) 


-3- 


-4- 


level  line  of  u;  i.e 


u(x  (a,  t ) ,  t)  »  a,  0  <  a  <  u  and  t  _>  T (1.17) 
where  T(a)  >_  0  18  the  first  tine  that  either  u  (t(a))  =  a  or  U(r(a))  =>  a,  and 


uQ(t) 


cAV(t)  *  <u  (t)  -  U(t)  )  3t  <ot't>da  ' 


(1.18) 


U(t) 


then  we  show  that 


lim  c„„(t)  »  c 

t-  AV 


Uf(u_) 

u 


(1.19) 


Equation  (1.19)  is  a  weak  statement  about  the  convergence  of  the  solution  of  the  Rieraann 

Problem  to  a  traveling  wave;  it  states  that  the  average  speed  of  propagation  of  the  level 

lines  of  u  converge  to  the  speed  at  which  the  traveling  wave  propagates, 
uf (u_) 

When  c  =*  —  >>  f'(0)  our  asymptotic  results  are  of  a  different  character.  Here 

we  show  that 

u(t  -  x/c)  <  u  <  u_,  cu(t  -  x/c)  <  o(x,t)  <  yf(u_),  and 

f(u(t  -  x/c))  -  cult  -  x/c)  <  ${x,t)  <  (1  -  y)f(u_)  (1.20) 

on  0  <  x  <  ct  and  t  >  0  and  that  the  shock  x  ~  s(t)  satisfies 


ct  <  s(t)  <  ct  + 


-J _ f 


-  (u 


u)(f ’ (u)  -  c) 


c ( c  -  f • (0) ) 


u(c  -  - 


!»f  (U)  1 


du 


(1.21  ) 


The  function  u(*)  is  the  u  component  of  the  traveling  wave  defined  earlier  in  item  (b) 
on  traveling  wave  solutions.  It  satisfies  (1.4)  on  T  >  0  and  the  initial  condition 
(1.7).  The  inequalities  (1.20)  and  (1.21)  also  guarantee  that 

f  (u  ,uf(u  ),(1  -  w) f (u  )),  X  <  c 

lim  <u,0,$)<Xt,t)  ~  (1.22) 

t**  |  (0,0,0) ,  X  >  c  . 


1 
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in 


2.  A-Priori  Estimates 
A.  General  Remarks 

In  this  section  we  develop  a-priori  bounds  which  guarantee  that  the  upper  bound  for 
the  interval  of  existence,  Tmax,  is  plus  infinity. 

We  confine  our  attention  to  the  characteristic  equations 


and 


<>t  +  f  ’  (u)0x  +  (1  -  u)o  -  u$  »  o  , 
$  -  (1  -  u)o  +  m4>  =  0  , 


(o) 

(♦> 


(o  +  c)  <==>  <y  +  #  <=  f(u) 


(u) 


which  we  insist  hold  in  0  <  x  <  s(t)  and  0  <  t  <  Tmax>  a  also  satisfies  the  boundary 
condition : 


-6- 


and  0  <  X  <  a ( T  ) 
—  max 


(2.2) 


and 

>  0,  0  <0,  and  $  <  0,  T(x)  <  t  <  T 

t  x  x  —  —  max 

The  function  t  -  T(x)  is  the  inverse  of  x  »  s(t!  and  satisfies 


dT  f'1  (Sj_T(jxm 
dx  “  E(T(x)) 


and  T(0)  =  0  . 


O 


def 

Proof.  We  start  with  the  observation  that  (p,q)  »  *°t'^t'  satisfies 


Pt  +  f'(u)px  -  k(u)(p  ¥  q)^  ¥  ( I  -  li)p  -  pq  «  0  , 


qt  -  <1  -  p)p  +  pq  =  0 


0  <  x  <  s(t)  , 
0  <  t  <  T„-„ 


p(0,t)  »  0  and  lim  q(x,t)  «  (1  -  p)E(t)  >0,  0  <  t  <  Tm 

x-»s(t) 
x<s(t) 


,-1 


Again,  the  u  «  f  (q  $)  and  u  *  k(u)  is  definea  by 

k(u)  d2f  f-(u)/(f >{u))2  . 


Lemma  2.1.  Suppose  p(x,t)  >  0  for  T(x)  <  t  <  Tmax  and  0  <  x  <  s(Tmax).  Then 
q(x,t)  >  O  for  T(x)  <_  t  Tfflax  and  0<_x<  S(Tmax).  □ 

Proof.  The  lemma  follows  from  the  identity 


q(x,t)  *  (1  -  H  )£  ( T(x))e  V  (  T(x"  * 


(1  -  M)  /  e'W(t  3,p(x,s)ds 
T(x  ) 


and  the  hypotheses  on  p  and  X .  □ 

Lemma  2.2.  p(x,t)  >  0  fc  T(x)  <  t  <  T  and  0  <  x  <  s(T_  ).  □ 

Proof.  We  start  with  the  observation  that 

lim  (pt  +  f 1  (u)px)(x,t)  >  k(u0(t))(l  -  u)2S2(0)e"2ut  +  uU  -  u ) K C 0 ) e “ u fc 

*  k(uQ(t))(1  -  ii)2y2f2(u  )e  2Wt  ¥  p2(  l  -  u)f(u  )e  Ut  > 
and  thus  we  are  guaranteed  that  p  is  positive  in  some  domain  T(x)  <_  t  <_  Tfflax  and 
0  <  x  <  r.  We  now  assume  the  lemma  i3  false.  Then  there  is  a  first  time,  t«,  and 
point,  x«,  with  T(x.)  <_  t„  _<  TBax  and  0  <  X,  <  s(Tmax)  such  that  p(x.,t,)  -  0 
(pt  ¥  f '  (u >px I  (x,, t„)  <_  0.  Rut,  (2.4),  implies  that 


(2.3) 

(2.4) 

(2.5) 

(2.6) 


0 

and 
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(Pt  +  f' (u)px)(x,,t,)  -  k<u(x,,t,))q2(x,,t(,)  +  uqtx^.t^) 
and  thus  Lemma  2.1  implies  that  the  right  hand  side  of  the  la3t  equation  is  positive.  This 
is  the  desired  contradiction  and  establishes  the  lemma.  O 

The  identities 

u  =  f  1 ( o  +  $)  and  f'(u)ut  *  p  +  q  (2.7) 

together  with 

f'(u)  >0,  p  >  0,  and  q  >  0  (2.8) 

imply  that 

-ut  »  a  <  0,  T(x)  <  t  <  T  and  0  <  x  <  s(T  v)  .  (2.9) 

t  x  —  —  max  —  max 

If  we  now  exploit  the  Identities 

♦(x,t)  •  (1  -  \i)  j  e_P*t_S'o(x,s)ds  , 

T(x ) 

*  «  -(1  -  u)T'(x)S(T(x))e*W(t‘T(x))  +  (1  -  v)  \  e_1J(t'8,o  (x,s)ds  , 

**  __ .  .  X 

T(x) 

and 

T'(x)  a  f  1  (I  (T(x) )  )/E  (T(x) )  and  S  (T(x) )  >  0,  0  _£  X  <  Tmax 

and  (2.9)  we  obtain 

<t>x  <  0,  T(x)  £_  t  £  Tmax  and  0  £  x  <  s(Tmax)  .  (2.10) 

This  concludes  the  proof  of  Theorem  2.1.  □ 

An  immediate  corollary  to  theorem  2.1  is 
Theorem  2.2.  Under  the  conditions  of  Theorem  2.1  the  following  bounds  prevail 

0  <  o  £  uf(u  ),  0  £  t  <  (1  -  u)f(u  ),  and  0  <  u  <  u_  (2.11) 

for  0  <  x  <  s(t)  and  0  <  t  <  T_„„.  □ 

3.  The  case  where  u_  is  small. 

We  now  focus  on  obtaining  upper  bounds  for  (o  ,$  )  «  (p,q)  and  lower  bounds  for 

(dx»$  )•  These  estimates  will  rely  on  the  assumption  that  u_  is  small.  For  simplicity 

we  shall  assume  that  's(u)  f " (u ) /( f '  ( u ) ) 2  satisfies 

Ik 

^  <u)  <  0  .  (2.12) 

-8- 
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This  assumption  guarantees  that 


kQ  =  k<0)  _>  k(u),  0  <_  u  <_  u_  . 

The  assumption  (2.12)  is  not  essential;  if  we  abandon  it,  the  same  estimates  obtain 

provided  we  redefine  kQ  by  kQ  max  k(u). 

0<u<u 


The  inequalities  0  uf(u  )  and  0  <  $  and  the  equation  q  “  ♦  »  (1  -  M)0 

imply  that 

q  £  M(1  -  u)f(u  ),  0  _<_  x  _<  s(t)  and  0  _<_  t  <_  Tmax  . 

Equations  (2.4)^  and  (2.5)  in  turn  imply  that 

pt  +  f  ■  <u)px  _<_  kQ(p  +  u( i  -  u)f(u_))2  -  (i  -  u)p  +  m2(i  -  u>ft u_>, 

0  <  x  <  s(t)  and  0  <  t  <  T 

InaX 

and 


p(0,t)  >  0,  0  <  t  <  T 


Moreover,  if 


flu  )  < 


(1  -  M) 


-  4uk 


then  (2.14)  may  be  rewritten  as 

pt  +  * '  <u)px  ±  k0<P  “  P_MP  *  P+).  0  <  X  <  s(t)  and  0  <  t  <  Tnax 

-/he  re 

0  <  p_  «  ( 1  }  -  1/(1  -  U)f(u  )  -  /(I  -  u)2  -  4kri)i(  1  -  u)f(u_)  < 


2k  „ 


2k, 


(1  -  li) 
2k„ 


U  (  1  -  |l)f  (u 


2k 


471 


(1  -  V) 


-  4kQu(1  -  u)f(u_) 


hn  immediate  consequence  of  (2.15)  and  (2.17)  and  the  bound  q  _<  u(1  -  y)f(u  )  is 
Theorem  2.3.  If  u_  satisfies  (2.16),  then 

A  <  y(1-u)f<u  )  and  a.  <  -  w(l-u)f(u  )  -  477m)2  -  4k  u(1-u)f(u  )  . 

t  -  -  t  -  2K0  ZK0 


Moreover,  the  identities 


and 


<1*  =  -ut  =  -(p  +  q)/f'(u) 
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t 

t 

t 


L 


(2.1  3) 


-  ut 


(2.14) 


(2.15) 


(2.16) 


(2.17) 


(2.18) 


(2.19) 


(2.20) 


»x  -  -(  1  -  g)f"1(S(T(x)))e‘ll(t'T(x,: 


-u( t-s  )„ 


{ 1  -  U)  /  e"M'"'°'<Mx,8)d8 


(2.21  ] 


T(  x ) 


imply  that 


-  Wf' To7  ^ _( 1  -  u)  ♦  -  11)2  -  4V°  "  u)f(u_>  )  (2.22) 


and 


♦  >_  - ( 1  — U  )u  e 


-U( t-T(x) ) 


■■■■I-lTm  (1_u  ‘  -  4k  u(i-w)f(u_)  )(i-e‘u<t‘T(x)))  . 


(2.23) 


2kotif'(0) 


d3f 


Remark  1.  If  we  assume  that  f  satisfies  the  additional  hypothesis  — -  <  0,  then  we 

du3 


find  that  if  (2.16)  holds  the  following  inequality  must  also  prevail: 


Uf(u  ) 


<  f’(0) 


(2.24) 


d3f 


Thus,  if  (2.16)  holds  and  — -  <_  0  the  traveling  wave  associated  with  the  states 

du3 


(u  ,  uf ( u  ))  and  (0,0)  at  x  equal  minus  and  plus  infinity  respectively  is  the  smooth 
one  described  in  Section  1.  C 

The  results  of  Theorems  2. 1-2. 3  imply  that  if  (2.16)  holds,  then  the  functions 


o ,  ♦,  o  »  $  ,  o  ,  end  $  are  uniformly  bounded  independently  of  T  so  long  as 
t  t  x  x 


£(t)  >  0  for  0  <  t  <  T  ,  Thus,  if  (2.16)  holds,  the  only  way  for  T„,„  to  be  finite 

—  max  DlaX 


if  1  satisfies 


lim  1( t)  =  0  . 


(2.25) 


t*T 


def 


Theorem  2.4.  The  function  Eft)  =  lim  d(x,t)  cannot  satisfy  (2.25)  for  any  finite 

x>s(t) 


T  .  D 
max 


x<s( t ) 

Proof.  Suppose  the  Theorem  is  false.  Then,  there  is  a  first  finite  time,  T  >  0,  such 
that  E { t )  >0,  0  _£  t  <  T  ,  and  (2.25)  holds.  To  show  this  cannot  happen  we  observe 
that 


—  ■~:-:WilhSYhr'*‘h  i  V*  '"i  .  v:  » 
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dl 

dt 


-  f ' (f-1 (£ ) ) ) 


(2.26) 


■  »  (Ot  +  S0X )(»(t),t) 


-(1 


-  ll)I(t)  +  0  ( S ( t ) 
X 


,t) 


t  Ut) 
f-’dl 


where  o  (s(t),t)  =  lim  o  (x.t).  The  a-priori  bound 
x*s(t)  * 
x<s(t) 

-  jJ— fV('oT  l  1  "  **  "  4l  -  u)2  -  4k0u(1  "  u)f(u  )  )  <_  °x  <  0,  0  _<  x  _<  s(t)  and 
0  t  <_  Tmax,  together  with  the  uniqueness  theorem  for  ordinary  differential  equations, 

applied  to  (2.26),  imply  that  if  lim  £(t)  *  0,  then  I(t)  S  0  for  0  <  t  <  T 

tr*T~ 

max 

This  contradicts  the  fact  that  lim  E(t)  «■  pf(u  )  >  0  and  establishes  the  Theorem.  □ 

t+0* 

Summarizing  the  results  of  this  subsection  we  obtain 

Theorem  2.5.  If  u  satisfies  (2.16),  then  T  »  +»,  and  the  bounds  of  Theorems  2.1- 

2.3  obtain.  G 

C.  The  case  where  u_  is  large. 

He  now  turn  to  the  case  where  u_  is  large.  We  start  with  a  lower  bound  estimate  for 
the  shock,  x  »  s(t).  This  estimate  is  valid  for  all  u_  >  0,  but  it  is  only  useful 
when  u_  is  large. 

Theorem  2.6.  The  shock  curve  x  =  s(t),  0  _<_  t  t_  Tmax,  satisfies 

uf(u  ) 

s(t)  >  ct  where  again  c  -  — ~  .  D  (2.27) 


Proof.  The  conservation  law  ut  +■  =  0,  boundary  condition  o(0,t)  =  uf(u  ),  and 

Rankine-Hugoniot  condition  (RH)  imply  that 

.  s(t)  s(t) 

—  /  u(x,t)dx  =  uf(u  )  and  /  u(x,t)dx  =  uf(u  )t  .  (2.28) 

dt  0  0 

The  upper  bound  u(x,t)  <  u  for  0  <  x  <  s(t)  and  0  <  t  <  T_.„  combines  with  (2.28)- 
to  yield  the  desired  result*  i  i 
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In  the  remainder  of  this  subsection  we  shall  assume  that 


Remark  3.  The  triple  (£,£,u)  defined  by  (2.33)  solves  (o),  (♦),  and  (a)  and  satisfies 
the  initial-boundary  conditions 

(£,£)(x,0)  -  (cu,f(u)  -  cuX-x/c),  0  <  x  <  cT(c)  (2.37) 


and 

£(0,t)  m  cult)  and  £(c(t  +  T(c)),t)  ■  f(u_(c))  -  c£_(c)  <0,  t  >  0  .  (2.38) 

Moreover,  for  any  curve  x  »  p(t)  satisfying  ct  <  p(t)  <  c(t  +  T(c))  we  have 

£(p( t ), t)  <  0. 

Proof  of  Theorem  2.7.  We  first  observe  that  if  (g,$,u)  solve  (o),  ($),  (u),  (BC),  and 
(RH),  then  the  properties  of  u(«)  and  the  fact  that  ct  <  s(t)  imply 

£( 0, t)  -  cult)  <  0(0, t)  -  uf(u_),  0  £  t  <  Tmax  ,  (2.39) 

and 


£(ct,t)  -  0  <  *(ct,t),  0  £  t  <  Tmax  .  (2.40) 

Lemma  2. 3.  Suppose  £<o  on  0  £_  x  <  ct  and  0  _£  t  <  Tmax.  Then  £<$  on  0  £_  x  £  ct 
and  0  £  t  ±  Tmax.  O 

Proof.  The  lemma  follows  from  (2.40)  and  the  identities 


$ (x, t)  »  ♦(x,x/c)e-lJ*t~1C|/c'  +  (1  -  u)  /  e"u(t"s^a(x,s)ds 


x/c 

t 


>  J(x,x/c)e-li^  t_x/c'  +  (1  -  y)  /  e~P(  t-s  )(j(x,s)d3  '  (2.41) 

x/c 


and 


t 

£(x,t)  -  (1  -  u)  /  e'U'  ”B  £(x,s)ds  .  □  (2.42) 

x/c 

The  preceding  lemma  together  with  (2.39)  imply 

£(0,t)  <  *(0,t),  t  >  0  (2.43) 

and  (2.39),  (2.40),  and  (2.43)  are  sufficient  to  guarantee  the  Theorem  is  true  in  some 
neighborhood  —  <  t  <  T_„„  and  0  <  x  <  e  with  the  caveat  that  $  »  $  «  0  at  the 

C  “  i’ll*  X  *•  —  — 

origin.  We  now  suppose  the  Theorem  is  false.  Then  there  is  a  first  time  t,  and  point 
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x.  with  —  <  t  <  T  ami  e  <  x*  <  cT  such  that  -  <7(xl(t,).  Moreover, 

*  c  —  — -  ma  x  —  ilia  x  #  •  —  *  * 

at  (x,,t4)  the  following  inequalities  hold: 

3  3 

+  f '  (u)  j^)(o-o)(xt,  t^)  £  0,  (^-^Hx^.t^)  >  0,  and  (u-u_)  (x,,  t„ )  >  0  .  (2.44) 

The  latter  two  inequalities  follow  from  Lemma  2.3  and  the  relations  u  =  f  1  (o  $)  and 

-  f_1(£  ♦  _£)•  4  little  algebra  shows  that 

+  f'(u)  |j)|o^)(x,,t,)  -  u($-*)(xt,t#)  +  ^(Xj.tJlfMa)  -  f 1  (u)  Kx^.t,  ) 

and  this  identity,  together  with  (2.44),  the  hypothesis  f"  >  0,  and  the  fact  that 

o^(x4,t4)  »  -c  jj-i  (t  ■  t,  -  x/c)  <  0  yield  the  contradictory  Inequality 
3  3 

+  f  •  (u)  5^0  (°  -  >  ®*  This  completes  the  proof  of  the  Theorem.  □ 

Our  next  task  is  to  obtain  an  upper  bound  for  the  shock  curve  x  =*  s(t).  We  start 
with  the  identity 


s(t) 

/  u(x,t)dx 
0 


Uf(u  )t 


ct 

/  u_dx 
0 


(2.45) 


which  we  rewrite  as 


8(t)  ct 

j  u(x, t)dx  «  /  (u_  -  u(x,t))dx 

ct  0 


(2.46) 


If  we  now  exploit  the  inequality  uy  <  0  and  (2. 30) j  we  obtain 


s(t)  ct  .  t 

u(s(t),t)(s(t)-ct)  <  /  u(x,t)dx  =,  /  (u_-u(x,t) )dx  <—  /  (u_-u(x))dt 


ct  0 

u ( T )  (u  -u)(f(u)-c) 


if 


ir 


u(c-uf (u)/u) 


t  -  (u_-u)(f'(u)-c) 
du  K  c  f  ~~ u(c-uf(u)/u)  du 


(2.47) 


v^here 


*<uj  Mf(u  ) 

u(s(t),t)  ■  lira  u(x,t)  «  f  (-(t))  and  c  - — — —  •  — - - 

x+s(t)  * 

x<s(t) 


(2.40) 


and  u(«)  is  the  function  defined  in  (2.3U.  The  Rank ine-Hugon lot  equation 
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\  ~i 

dt 


£(u(a(t),t) ) 
u(s(t ), t) 


together  with  (2.47)  imply  that  y(t)  »  e(t) 


ct  satisfies 


dy 

dt 


< 


-(c 


f"  -  (u  -  u)(f(u)  -  c) 

f.,o»  ♦  (-=5E/ 

u 


u(c  -  uf(u)/u)  du'  y 


(2.49) 


and  y(0)  =»  0.  The  number  f*  „  is  given  by  f"  „  -  max  f**(u).  An  immediate 

max  max 

0<u<u 

consequence  of  (2.49)  and  y(0)  ■  0  is 

Theorem  2.9.  The  shock  curve  x  »  s(t)  satisfies 

s(t)  <  ct  +  c5(c),  0  _<_  t  Tnax  ,  (2.50) 

where 


, ,  .  max  r 

He)  .  - -  ; 

2(c  -  f'(0))c  u4 


-  (u_  -  u)(f*(u)  -  c) 
u(c  -  uf(u)/u) 


du 


□ 


(2.51 ) 


Remark  4.  If  we  restrict  our  attention  to  the  case  where  f(u)  »  f'(0)u  ♦ 


V 


it  is 


easily  established  that  the  numbers  T(c)  of  (2.35)  and  4(c)  of  (2.52)  satisfy 

L)  log (  (  2-p )c-yf 1  ( 0)  1  p(c-f  ■  ( 0) ) 


_ 2|f.  U  (c- 

t ( c )  2(c'f  (0))c  ^  ‘  T  Tc-iTf  *To) 

7(0 - 


2<1-u)c 


2(c-uf»(0)) 


log  2) 


(C.f(0),loq(^miJ} 


>  1  (2.52) 


uniformly  on  f*(0)  <<  c  and  0  <  1  -  u  <<  1. 


□ 


Motivated  by  the  results  of  Remark  4  we  shall  make  the  following  additional  assumption 

T(c) 


5(c) 


>  1 


(2.53) 


uniformly  on  f'(0)  <<  c  and  0  <  1  -  u  <<  1.  Q 

Using  the  results  of  Remark  3  and  the  arguments  used  to  establish  Theorem  2.7,  one 
easily  obtains 


yf  (u  ) 

Theorem  2.9.  If  (2.53)  holds,  then  for  f'(0)  <<  c  -  — - -  the  conclusion  of  Theorem  2.7 

are  valid  on  the  domain  0  <  x  <  s(t)  and  0  <  t  <  T _ .  □ 

Our  final  task  in  this  subsection  Is  to  obtain  upper  bounds  for  p  »  ct  and  q  » 
and  lower  bounds  for  and  Throughout,  we  shall  assume  that 

k(u)  f " (u)/( f 1  (u)  )2  in  decreasing  on  u  >  0.  We  also  assume  that  (2.53)  holds  and 
that  o  >  f ' ( 0)  is  large  enough  so  that  the  conclusions  of  Theorem  2.9  hold.  Finally,  we 


-15- 


shall  assume  that 


satisfies 


.  <ief  uf'(u)  -  f(u)  .  _  _  . 

>(u)  '  - TUT) - >  0>  0  <  u 


|1  (u)  >  o,  o  <_*  . 


Our  goal  is  to  show  that  subject  to  the  above  hypotheses,  there  is  some  0  <  X  <  u 


such  that  the  functions  p  «  and  q  »  satisfy 


p(x,t)  «e  Mt-T(,t),P(x,t),  and 


q(x,t)  »  e”A(t_T<X> ^Qtx.t) 


T(x>  ±  *  <-  Tmax  ' 
0  i*  <.stTmax> 


where  P  >  0  and  p  >  0  are  bounded  independently  of  Tma){.  Again,  t  ■  Ttx)  is  the 
inverse  of  the  shock  x  «  stt). 

It  is  a  relatively  simple  matter  to  show  that  if  p  and  q  are  given  by  (2.56)  and 


(2.57)  and  satisfy  (2.4)  and  (2.5),  then  P  and  Q  satisfy 


Pt  ♦  f 1  (u)Px  _<_  k(u)e 


(P  +  Q)  -  (1  -  u  Xj (u) )P  ♦  uQ  ,  (2.58) 


Qt  -  (1  -  U)P  -  (U  -  X)Q  ,  (2.59) 

in  the  domain  0  <  x  <  s(t)  and  0  <  t  <  Tmax,  and 

p(0,t)  -  0  and  lim  Q(x,t)  -  (1  -  U)I(t),  0  <  t  <  Tmax  .  (2.60) 

x»s(t) 
x<s(t) 


The  results  of  Theorem  2.1  guarantee  that  P  and  Q  are  positive  on  0  <  x  <  s(t)  and 

0  <  t  <  T___  while  the  results  of  Theorem  2.9  and  (2.53)  guarantee  that 

maX 

u(x,t)  >  u(c)  >0,  0  <  x  <  s(t)  and  0  <  t  <  T___  (2.61) 

—  —  «■—  —  max 

where  again  u_(c)  is  defined  by 

f(u_(c))  -  c  (2.62) 

and  11m  u_(c)  -  ¥*,  The  Inequality  (2.61),  when  combined  with  the  hypotheses 

<  0  and  >  0  , 
du  du 

and  the  results  P  >  0  and  Q  >  0  for  0  <  x  <  s(t)  and  0  <  t  <  Tnax  imply  that  P 
and  Q  satisfy  the  Inequalities 
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(2. *3) 


pt  ♦  f'(u)Px  <  k_(p  +  g)2  -  O  -  u  ♦  ^i)P  ♦  up  , 

Ot  -  (1  -  u)P  -  (M  -  X)0  .  (2.64) 

on  the  sane  domain.  The  numbers  j<_  and  j_  are  defined  by 

k_  -  k(u_(c) )  and  0  <  j_  •  j  (ji(c) )  .  (2.65) 

Remark  5.  If  me  restrict  our  attention  to  the  case  f(u)  »  f'(0)u  ♦  kQu2/2,  then 


k(u) 


7.  j<u) 


V 


(f*<0)  ♦  kQu) 


( 2f 1 ( 0)  +  kQu) '  - 


u(c ) 


(c  -  f(0)) 


and  lim  (k(u(c ) ), j (u(c) ) )  ■  (0,1).  n 

c*» 

Motivated  by  the  last  remark  we  shall  assume  that  k(*)  and  j(‘>  satisfy 


lim  ( k < jj_( c ) ) ,  j  (u( c )  > )  -  (0,1)  . 

c  ►» 


(2.67) 


We  now  turn  to  the  inequalities  (2.63)  and  (2.64).  Throughout,  we  assume  that 
0  <  X  <  U.  It  is  easily  checked  that  the  set  P  _>.  0  and  0  _>_  0  such  that 

k_(P  +  o)2  -  (i  ♦  -  u)p  ♦  up  ■  0  is  given  by 

q  »  Q(p)  dSf  -p  ♦  /u2  *  4k ( t  ♦  \J.)  -  U,  0  1  p  <  <'  +  M  "  (2.66) 

and  that 

MP  ♦  Q)2  -  (1  +  -  U)P  ♦  uO  <  0,  0  £  0  <  P(P>  and  0  <  P  <  ( 1  +  -  M >  •  (2.66) 

It  is  also  a  simple  matter  to  verify  that  if  (u  -  X>i  -  (1  -  u)  >0,  then  the  right  hand 
sides  of  (2.63)  and  (2.64)  vanish  simultaneously  at  where 


Mu  -  A)((U  -  Mj_  -  (1  -  w>) 


eq 


(1 


-  X)^ 


>  0 


(2.70) 


and 


*eq 


X  ( 1  -M)((M-X)j_-(!  -  w  )  ) 
(1  -  \)2k 


>  0 


We  now  choose  X  so  as  to  maximize  p^^(X).  The  result  is 


Mi  -  U  -  R)  M(1  t  i>  -  1 


max  (2  -  »i)i  ♦  1  -  M  (2  -  M>i  ♦  1  -  M 
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(2.7)  ) 


(2.72) 


Substitution  of  (2.72)  into  (2.70)  and  (2.71)  yields 


p  (n,j)  dSf  P  (X  ) 
eq  «-  eq  max 


(Uj_+1  +M )  (irj_  -  ( 1  -U )  <  2<  2U-1  )  -  ( 1  -j_)  (  2U-1  *-U ( 1  +J_) ) ) 
4((2-U)i  +  (l-lOHl+i)2* 


(2.73) 


0  (u.i)  dSf  o  (X  ) 

eq  —  eq  max 


(Wi  -  <1-U))(2(2U-1)  -  ( 1  -j.)  (  2P-1  +M(  1  +j_) )  ) 
- - ,  (2.74) 

4(1 +i>  \ 


and  these  formulas,  together  with  (2.67),  imply  that  if  0  <  1  -  u  <<  1  i  «  c, 

then  P-f^.j.)  ond  C>  (M,j)  are  both  positive;  in  fact  we  have 

lim  ()Lp_n < w ,  j  ) , kO  ( U ,  j  ) )  »  (|,  -j)  .  (2.75) 

c  h»  eq  —  —  eq  —  88 

u»r 

Our  final  task  is  to  obtain  conditions  which  quarantee  that  P  <  p  <u,j)  and 

'tq  jL- 

0  <  0  (U,j>  for  0  <  x  <  s(t)  and  0  <  t  <  T  The  fact  th.it 

eq  max 

lim  o(x,t)  =  £  ( t )  yf(u  )  for  0  <  t  <  Tmax  together  with  the  boundary  condition 

x*s(t) 
x<s( t) 


lim  Q(x,t)  =  (1  -  u)£(t)  guarantee  that  Q  <  p  <M,1)  on  x  =»  s(t)  provided 
x*s(t)  eq 

x<s(t) 

(Ui-(1-U))(2(2u-1)-(1-JL)(2u-1+lJ(l+i))) 


u('-u)f(u  )  <  p  (u,i)  or  u  ( 1-u  )kf  (u  )  < - 

-  eq  —  - 


-.  (2.76) 


4(1+i>‘ 


uf (u  ) 

Moreover,  the  hypotheses  on  f,  k,  and  j  quarantee  that  for  c  »  f'(u(c))  - - — 

—  u 

sufficiently  large  the  last  inequality  holds  provided  0  <  1  -  M  <<  1.  With  these 
preliminaries  we  are  now  in  a  position  to  prove 

Theorem  2.10.  Suppose  the  functions  k  >  0  and  j  >  0  satisfy  dk/du  <  (1  and 
dj/du  >  0  and  that  (2.67)  holds.  Suppose  further  that  f'(0)  <<  c  and  0  <  1  -  u  <<  1 


are  such  that  (i)  (2.53)  holds,  (ii)  that  the  functions  p  (u,j)  and  p  (u,j)  defined 

eq  eq  *- 

in  (2.73)  and  (2.74)  are  both  positive,  and  (iii)  that  (2.76)  holds.  Then  the  functions 

-  (1  -  H) 

P  and  p  defined  in  (2.56)  and  (2.57)  with  X  «  — - — - - - -  >  o  satisfy 

\2  -  U  f  Vi  -  U ) 
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(2.77) 


0  <  P(x,t)  <  and  0  <  Q  <  Q  (u,i> 

for  T(x)  <_  t  _<_  Tmax  and  0  <  x  <_  s(Tmax).  n 

Proof .  The  lower  bounds  follow  from  the  results  of  Theorem  2.1  and  as  previously 

observed  Q  satisfies  the  desired  upper  bound  on  the  boundary  x  «  s(t),  0  <  t  <  T  . 

—  —  max 

The  identity  lira  Q(x,t)  -  u ( 1  -  u)f(u  )e_{ U-X)t  together  with  ( 2. 76 )  and 
x+0 

2(1  -  U)j_ 

<M  -  1)  «  -■  —  -  >  0  imply  that  0(0,  t)  <  for  0  <_  t  <_  Taax.  The 

differential  inequality  (2.63),  together  with  (2.64)  and  the  results  summarized  in  (2.68) 

and  (2.69)  imply  that  the  upper  bounds  cannot  fail  in  T(x)  <  t  <  T  and 

—  —  max 

0  <  x_<_s(Tmax)  and  this  concludes  the  proof.  O 

Lower  bounds  for  and  9^  are  readily  obtainable  from  the  upper  bounds  for  o 

and  ♦  .  The  results  are  that  if  the  hypotheses  of  Theorem  2.10  hold,  then 


.  _ -X ( t-T(x) )  ‘V*1  ♦  W*” 

-*  rco) 


(2.78) 


and 


♦.  >  -(1-U)U  e-‘i(t-T(5t>)  -  TtExT  <peq(w'i,  ♦  0ea(M'i))(®"X<t‘T<X>>-«’li<t"T<X)))  (2-79) 


eq 


Mi  -  (1  -  M) 

for  T|xl<t<tI11(  and  0  <_  x  <_  s(Tnax)  with  1  -  jy— »  (t  .  >  0  and 

2(1  -  M>2 

(u  -  X,  .  ______  >  o.  Finally,  since  the  bounds  of  Theorems  2.7  throuqh  2.10 

and  the  inequalities  (2.78)  and  (2.79)  are  independent  of  #  we  obtain 

Theorem  2.11.  If  the  hypotheses  of  Theorem  2.10  hold,  then  T  «  +*>  and  the  bounds  of 
-  max 

Theorems  2.7  through  2.10  and  the  inequalities  (2.78)  and  (2.79)  hold  with  T  «  ♦*.  □ 

max 

This  concludes  Section  2. 
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3.  Asymptotic  Results 

A.  The  case  where  u_  is  small 

In  this  section  we  shall  assume  that  (2.12)  and  (2.13)  hold.  We  shall  also  assume 

that  u_  is  small  enough  so  that  (2.16)  is  valid.  Our  goal  is  an  asymptotic  estimate 

for  U(t)  lim  u(x.t).  The  Rankine-Hugoniot  conditions  ( 1 . 1 1 )-( 1 . 1 3)  imply  that  U 

x*s(t) 
x<s(t) 


satisfies 


(U) 


and  u(0)  =  Ug  , 


(3.1 ) 


where 


ut(s(t),t)  -  lim  u  (x,t)  -  (<1  -  W)  -  /(I  -  u)2  -  4k  u(1  -  |i)f(u  )  )  ,  (3.2) 

x*s(t)  "o  1  ’  0 

x<s(t) 


and  0  <  un  <  u_  is  the  unique  solution  of 


f(u0)  -  uf(u_) 


The  fact  that  kgf'(O)  ■  f"(0)/f'(0)  and  the  inequalities 
(1  -  M) 


-  /(I  -  U)2  -  4k  |i(1  -  u)f(u  )  <  4k  uf (u  >  <  (1  —  u ) 


0  -  0 
when  combined  with  (2.16),  (3.1),  and  (3.2)  yield 

dU . .f  f( 0)  (.  f(U)  i  f(U) 


The  fact  that 


satisfies 


h(U) 


def  f(0) 


(1 


f CO)  ^  f(U) 


2f " (0)  v  Uf'(U);  f'(U) 


(3.3) 


(3.4) 


au  .  .f  fioi  i ,  tiui  i  riu)i  .  ,, 

dt  —  '  ”  W  ^2f"(of  1  “  Uf  * (U)  “  f '  (uf  and  U  0  “  u0  •  (3.5) 


(3.6) 


h(0)  -  0  and  1^-  (0)  -  -3/4  <  0 

guarantees  that  if  Ug  and  hence  u_  are  small,  then  u(t)  *  0  a9  t  ♦  in  fact  for 

Uq  small  lira  e^fcu(t)  »  0  for  any  0  <  X  <  -3/4.  In  the  special  case  where 
t»“ 

f(U)  -  f'(0)U  +  k0tl^/2,  it  is  easily  verified  that 


h(U) 


(3f'(0)U  *  2k  QU  ) 
4( f 1  ( oT  ♦  'k'uj 


<  0 


(3.7) 


and  thus  no  additional  restrictions  on  Ug  are  required. 
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i 


We  conclude  this  subsection  by  showing  that  when  u_  is  small  (in  particular  small 
enough  so  that  U(t)  *  0  as  t  *  ")  the  solution  to  (B),  (C),  and  (1.2)  converges,  in  a 


weak  sense,  to  the  traveling  wave  described  in  Section  1.  Specifically  we  shall  prove 
Theorem  3,1.  The  average  speed  of  propagation  of  the  level  lines  of  u  converge  to  the 
speed  of  the  traveling  waves  that  is  the  function 


uo(t> 


CAV(t) 


UQ(t>  -  U(t) 


J 

U(t) 


3JL 

at 


(a,t)da 


(3.8) 


satisfies 


lim  c 

t*» 


AV 


(t) 


Uf(u_) 

u 


(3.9) 


The  function  Ug(t)  is  defined  in  (1.9).  For  numbers  0  <  a  <  u_,  the  curve 
x  ■  x(a,t),  t  _>  T(a),  is  a  level  line  of  u;  that  is  it  satisfies  u(x(a,t),t)  -  a. 

For  0  <  a  <  Uq(0),  r(a)  is  defined  by  U(t(a))  «  o  and  for  uQ{0)  £  a  <  u_,  t(a)  solves 
uQ(t(a))  »  a.  □ 

Proof.  Our  first  task  is  to  show  that 


uf(uj  -  f(U(t)) 
cAV(t)  “  uQ(t)  -  U(tj  ’ 

To  obtain  (3.10)  we  note  that  the  defining  relation  u(x(a,t),t)  »  a 

u  x  »  1  and  u_  +  xtu„  «  0.  If  we  now  define  E(a,t)  ■  o(x(a,t),t) 
x  a  t  t  x 

Z  (a,t)  »  a  x  and  this,  combined  with  u„  =■  -a  and  x,.  «  -u.x  , 
a  x  a  txcto 

conservation  law: 


(3.10) 

Implies  that 
,  we  see  that 
yields  the 


xt(a,t)  »  to(o,t),  U(t)  <  a  £  u0(t)  .  (3.11) 

Equation  (3.10)  now  follows  from  (3.8)  and  (3.11)  and  the  limit  relation,  (3.9),  is  a 

consequence  of  (3.10)  and  lim  (U(t),uQ(t))  «  (0,u_).  n 

t*» 


B.  The  case  where  u_  is  large 

To  obtain  the  asymptotic  result 


lim  (u,  J  ,<>)  (Xt,  t) 


(3.12) 


C  uf(u_) 

(u  ,uf(u  ),(1  -  u)f(u  )),  X  <  c  «  - 

-  - 

I  uf(u  ) 

!  (0,0,0),  X  >  c  - — 

V  II 


11  that  is  required  are  the  lower  bound  estimates  (2.30),  the  inequalities  (2.27),  (2.51) 
and  (2.52),  and  the  assumption  that  Tmax  «  +"  (which  we  know  is  true  if  the  hypotheses  of 
Theorem  2.10  hold).  When  X  >  c,  the  curve  x  »  Xt  satisfies 

s(t)  <  Xt  (3.13) 


for  times  t  > 


2( X  -  c)(c  -  f '  ( 0 ) ) c 


-  (u  -  u) (f 1 (u)  -  c) 
u(e  -  uf(u)/u) 


du  T(X,c). 


Thus 


(u, J ,$) ( x, t )  2  0  for  s(t)  <  x  and  t  >  T(X,c)  establishes  (3.13)j. 
follows  from  (2.30)  and  the  fact  that  when  X  <  c 

lim  u(t  -  —  t)  ■  lim  u(T)  =  u  . 

t  *•»  C 


That  (3.1 3)1  holds 


(3.14) 


C.  Concluding  remarks  and  open  questions 

As  was  mentioned  earlier,  we  have  succeeded  in  showing  that  when  u_  is  small, 

cftv(t)  converqes  to  the  speed  of  the  traveling  wave.  One  would  also  like  to  know  if  this 

result  is  true  when  u_  is  large.  The  result  would  follow  from  (3.10)  provided  we  knew 

f 

that  when  u_  was  large,  IJ ( t )  =  lim  u(x,t)  converged  to  the  number  u«  (defined  by 

x^sft) 

x<s(t) 


f(u.)  ut(u_) 

— u —  *  — u - '  as  6  approached  infinity. 

is  equivalent  to  establishing  that  lim  s(t) 

t  w> 


At  the  moment  this  is  an  open  question  (and 


f(ut> 


Assuming  for  the  moment  that 


lim  n(  t)  »  y 
t>»  (  0  , 


uf(u  )/u  >  f (0)  ' 
uf (u  )/u  <  f'(0)  | 


and  hence  that 


yf (n_  ) 

lim  cfty(t)  =  — j -  we  wool  1  also  like  to  know  whether  o(u,t)  =  x^iu.t)  =  E^(q,t)  ( seo 


» 


for 


(3.11))  satisfies  limc(a,t) 

t+~ 


Uf(u_) 

a 


j  ut  <  a  <  u  when  uf(u  )/u  >  f'(0) 
|  0  <  a  <  u  when  yf(u  )/u  <  f 1 ( 0 ) 


Finally,  we  end  with  a  tantalizing  calculation  which  would  be  sufficient  to  guarantee 
Uf (u  ) 

that  lim  U(t)  =  u,  when - —  >  f'(0).  The  differential  equation 

(-M0 


dU 

dt 


ut<s(t),t)(l 


f  (U) 
uf  ~(u) 


together  with  the  identity 


(1 


M> 


f(U) 
f  •  <  U) 


c(U(t),t) 


f '<U(t))u  (s(t),t) 

U  (s(t), t)  +  (1  -  U)U(t) 


implies  that  T)  satisfies 


do 

dt 


(1  -  u)(c(U(t),t)0(t)  -  f  (IJ(t) ) ) 
(f’(U(t))  -  c(U(t),t)) 


and  0(0)  =  un 


(3.15) 


where  again  c(a,t)  =  xt(a,t)  »  t)  and  uQ  is  defined  in  (3.3).  tf  we  now  knew 

uf(u  )  -  f(0) 

that  cAV(t)  =  ~o(-t)  -  0(tj  satlsfied 

cRy(t)  <  c(0(t),t) 

then  the  inequality  uQ(t)  <  u_  and  (3.15)  would  imply  that  0  okeys  the  inequality 


do 

dt 


> 


(1  -  u) (uf (u_)u  -  f (0)U  ) 

Tf '  <U)  (u  •  0)  -  (U (u  )  -  f(0)  )  ) 


and 


0(0) 


(3.16) 


An  immediate  consequence  of  (3.16)  is  that  'J(t)  >  u„  for  all  t  >  0.  At  this  point  we 

would  be  able  to  conclude  that  lim  Oft)  =  u*.  The  result  would  follow  from  the  upper  and 

t  e* 

lower  bounds  for  the  shock  curve  (see  (2.50))  and  from  the  fact  that  — —  is  uniformly 

dt 

bounded  on  [0,*). 
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